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'''The  nonlinear  evolution  of  the  Free-Electron  Laser/Ubitron  amplifier  is  investigated  numerically 
for  a  configuration  consisting  of  helical  wiggler  and  an  axial  guide  magnetic  fields.  A  set  of  coupled  non¬ 
linear  differential  equations  is  derived  in  three-dimensions  which  governs  the  self-consistent  evolution  of 
either  the  TE  or  TM  modes  in  a  loss-free  cylindrical  waveguide  and  the  trajectories  of  an  ensemble  of 
electrons.  The  initial  conditions  are  chosen  to  model  the  adiabatic  injection  of  a  cold,  cylindrical^  sym¬ 
metric  electron  beam  into  an  interaction  region  in  which  the  wiggler  amplitude  rises  slowly  from  zero  to  a 
constant  level  in  ten  wiggler  periods  Both  self-field  and  space-charge  effects  have  been  neglected  in  the 
formulation,  and  the  analysis  is  valid  for  the  High-Gain  Compton  regime  of  operation.  Numerical  simu¬ 
lations  are  conducted  to  model  an  amplifier  operating  in  the  neighborhood  of  35  GHz,  and  for  electron 
beam  energies  of  250  KeV  and  1  MeV.  The  growth  rate  in  the  linear  regime  prior  to  saturation  is  found  - 
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to  be  in  substantia)  agreement  with  the  predictions  based  on  a  linear  theory  of  the  instability,  and  the 
saturation  efficiency  is  consistent  with  that  expected  on  the  basis  of  simple,  heuristic  phase-trapping  argu¬ 
ments  Substantial  enhancements  in  the  efficiency  are  found  to  occur  due  to  the  presence  of  the  axial 
guide  field 
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NONLINEAR  ANALYSIS  OF  UBITRON/FREE-ELECTRON  LASER 
AMPLIFIERS  IN  THREE-DIMENSIONS 


I.  INTRODUCTION 

The  Ubitron/Free-Electron  Laser  (FEL)  has  been  successfully  demonstrated  for  operation  over  a 
broad  frequency  range  from  the  microwave  to  the  optical  parts  of  the  spectrum,  and  shows  promise  as  a 
high  power  source  of  radiation  at  millimeter  and  submillimeter  wavelengths.1-11  The  distinction 
between  the  Ubitron  and  the  FEL  is  not  well-defined  in  the  literature.  In  our  work,  we  find  it  con¬ 
venient  to  distinguish  between  the  FEL  and  the  Ubitron  primarily  on  the  basis  of  the  electron  beam 
energy,  and  we  shall  refer  to  such  devices  as  Ubitrons  when  the  electron  beam  energy  is  less  than  or  of 
the  order  of  500  keV.  While  this  choice  of  energy  is  arbitrary,  it  is  motivated  by  the  fact  that  most 
devices  that  operate  at  energies  below  500  keV  are  characterized  by  operation  frequencies  that  occur 
close  to  the  waveguide  cutoff.  The  Ubitron12  therefore,  may  be  though  of  as  a  weakly  relativistic  FEL 
operated  as  a  microwave  tube,  and  we  emphasize  that  the  physical  mechanisms  for  the  FEL  and  the 
Ubitron  are  identical. 

The  motivation  for  the  present  work  is  to  develop  a  fully  three-dimensional  nonlinear  analysis  and 
simulation  code  of  the  FEL/Ubitron.  Theoretical  investigations  of  the  Ubitron/FEL  have  been  con¬ 
fined,  principally,  to  the  linear  regime;13-24  however,  a  fully  nonlinear  treatment  is  required  to  describe 
the  interaction  through  the  linear  stages  to  saturation.  We  judged  it  important  to  include  the  effect  of 
an  axial  guide  magnetic  field  on  the  interaction  in  the  analysis  because  (1)  axial  guide  fields  are 
included  in  many  experimental  configurations  to  confine  high  current  electron  beams  against  the  effects 
of  self-fields,  (2)  linear  analyses  of  the  interaction  have  shown  that  substantial  enhancements  in  the 
gain  are  possible  due  to  the  presence  of  the  guide  field,  and  (3)  nonlinear  simulation  in  one-dimension 
Manuscript  approved  July  22,  1985. 


1 


has  shown  corresponding  enhancements  in  the  saturation  efficiency.  Three-dimensional  effects  become 


i  important  when  the  transverse  quiver  velocity  due  to  the  wiggler  magnetic  Field  becomes  large.  In  this 

I  regime,  the  electron  displacement  from  the  axis  of  symmetry  also  becomes  large,  and  a  one- 

t 

dimensional  approximation  for  the  wiggler  magnetic  field  breaks  down.  In  addition,  the  transverse 
mode  structure  of  the  radiation  field  is  also  important  in  the  description  of  the  overlap  of  the  radiation 
I  field  and  the  electron  beam,  and  is  a  crucial  feature  of  the  analysis  when  the  interaction  occurs  in  the 

vicinity  of  the  waveguide  cutoff  in  Ubitrons. 

In  this  paper,  we  derive  a  fully  three-dimensional  nonlinear  analysis  of  the  FEL/Ubitron  for  a 
configuration  in  which  an  energetic  electron  beam  is  propagating  through  a  loss-free  cylindrical 
waveguide  in  the  presence  of  a  helically  symmetric  wiggler  and  a  uniform  axial  guide  field.  To  this 
end,  a  set  of  coupled  nonlinear  differential  equations  is  derived  which  self-consistently  describes  the 
evolution  of  both  an  ensemble  of  electrons  and  the  electromagnetic  fields.  Space-charge  fields  are 
neglected  in  the  analysis;  therefore,  the  treatment  is  applicable  to  the  High-Gain  Compton  (or  strong- 
pump)  regime.  The  nonlinear  current  which  mediates  the  interaction  is  computed  from  the  microscopic 
behavior  of  an  ensemble  of  electrons  by  means  of  an  average  of  the  electron  phases  relative  to  the  pon- 
deromotive  wave  formed  by  the  beating  of  the  radiation  and  wiggler  fields.  These  equations  are  solved 
for  the  case  in  which  a  monoenergetic  electron  beam  of  arbitrary  initial  cross-section  is  adiabatically 
injected  into  the  interaction  region.  The  adiabatic  injection  is  modelled  by  allowing  the  wiggler  field 
amplitude  to  increase  slowly  from  zero  in  ten  wiggler  periods.  The  finite  waveguide  geometry  is 
included  in  the  analysis  by  the  introduction  of  the  boundary  conditions  appropriate  for  either  the  TE  or 
TM  modes  in  a  loss-free  cylindrical  waveguide.  Thus,  all  transverse  and  finite  geometry  effects  are 
included  in  the  static  wiggler  and  radiation  fields  in  a  self-consistent  manner.  In  addition,  since  the 
problem  of  interest  is  that  of  a  Ubitron/FEI  amplifier,  only  single  wave-mode  propagation  is  con¬ 
sidered.  This  permits  an  average  over  a  wave  period  to  be  performed  which  eliminates  the  fast  time 
scale  phenomena  from  the  formulation.  The  resulting  equations  are  equivalent  to  a  kinetic  description 
of  the  problem,  and  result  in  a  great  increase  in  computational  efficiency  over  a  full-scale  particle-in-cell 


simulation. 


The  organization  of  the  paper  is  as  follows.  The  general  equations  are  derived  in  Sec.  II.  The 
numerical  solution  of  the  coupled  particle  and  Held  equations  is  given  in  Sec.  Ill  for  parameters 


appropriate  to  both  Ubitrons  and  FEL’s.  A  summary  and  discussion  is  given  in  Sec.  IV. 

II.  GENERAL  EQUATIONS 


The  physical  configuration  we  employ  includes  a  uniform  axial  guide  field  and  a  helically  sym¬ 
metric  wiggler  field  generated  by  a  bifilar  helix,  so  that  the  static  magnetic  Field  can  be  written  in  the 
form 


B(x)  —  B0ez  +  B„(x),  (1) 

where 


Bh,(x)  -  2Bw(z)[l[  (X)  e,  cos*  -  -|-/|(X)e#  sin*  +  /, (A )  ez  sin*],  (2) 

represents  the  wiggler  field  in  cylindrical  coordinates.  In  Eq.  (2),  Bw(z)  describes  the  wiggler  ampli¬ 
tude,  X  =  k„r,  \  =  9  -  kwz ,  kw  (=  2 ir/kw,  where  Xw  is  the  wiggler  period)  is  the  wiggler 
wavenumber,  and  /„  (X )  and  l„(k)  represent  the  modified  Bessel  function  of  order  n  and  its  derivative 
respectively.  The  adiabatic  injection  of  the  electron  beam  is  described  by  allowing  the  wiggler  ampli¬ 
tude  to  vary  slowly  in  z,  and  is  valid  representation  as  long  as 

-7-  lnfi„  «  kw. 
dz 

In  practice,  we  shall  allow  Bw(z)  to  vary  only  over  0  <  z  <  10xw,  after  which  it  shall  be  held  constant. 
In  the  simulation  it  is  assumed  that 


B,(z) 


Bw  sin2  (kwz/ 40)  ;0  ^  z  <  10XW 

Bw  ;z  >  10X* 


(3) 


Since  the  space-charge  fields  are  neglected,  the  boundary  conditions  at  the  waveguide  wall  may  be 
satisfied  by  expanding  the  vector  potential  in  terms  of  the  orthogonal  basis  functions  of  the  empty 
guide.  Thus,  we  write  the  vector  potential  of  the  radiation  field  in  the  form 


8/4  (x,  /)  “  X  T —  ,r)  e,  sina,  +  Jt(.kinr)  es  cos  a, 

/»o  */» r 


for  the  T£  modes,  and 


00  ,  /  kfa 

8A(x,  /)  -  £  Jt(klnr)  e,  cosa(  — - —  Ji(k,„r)  et  sin  a,  +  —r  Ji(khr)  et  sin  a,  ,  (5) 

;-0  "to r  * 


for  the  TM  modes,  where  for  frequency  to  and  wavenumber  k(z) 

a  1  =  f  dz'k(z)  +  19  —  ot  1.  (6) 

In  Eqs.  (4)-(6),  J,  and  J[  represent  the  regular  Bessel  function  of  the  first  kind  of  order  /  and  its 
derivative,  and  k,„  describes  the  cutoff  wavenumber.  In  the  case  of  the  T£,„  mode  kln  =  x'j Rt ,  where 
J[(x[n)  -  0,  and  is  the  waveguide  radius.  For  the  TMtn  mode  k,„  =  xtJRg,  where  J, (x,„ )  -  0. 
Thus,  there  is  an  implicit  limitation  on  the  magnitude  of  the  beam  currents  which  can  be  treated  self- 
consistently,  and  it  is  assumed  that  the  mode  amplitude  6/4/,,  (z)  and  wavenumber  k(z)  are  both 


slowly-varying  function  of  z  such  that  both 


In  8/4/*  (z)  «  k. 


In  k(z)  «  k. 


The  microscopic  source  current  can  be  written  as  the  following  sum  over  identical  particle  trajec¬ 


tories: 


sii  ,  L  ~  .  ,  8ft  t j  (z ,  Xf0 ,  yi0 ,  tm )  1 

t)  *  €flb  XfQ ,  I  /  \(  ’  ^ 

>«•  j  » vzj  iz ,  x^  ,  j'jio »  '  1 

where  L  is  the  length  of  the  interaction  region,  NT  is  the  total  number  of  electrons,  nb  is  the  average 
electron  density,  v*(z;  x*,>to,  ti0)  is  the  velocity  of  the  / th  electron  at  position  z  which  entered  the 
interaction  region  (i.e.,  crossed  the  z  **  0  plane)  at  time  tio  and  transverse  position  (x/0,  >>,<,),  and 


t  1  (z ,  xio  ,  ylo ,  (l0 )  *  tl0  H 


vri(z  ,  xto ,  y,0 ,  tl0 ) 


The  system  is  assumed  to  be  quasistatic  in  the  sense  that  particles  which  enter  the  interaction  region  at 
times  t0  separated  by  integral  multiples  of  a  wave  period  will  execute  identical  trajectories.  As  a  result. 


v,(z;  Xj0 ,  yi0 ,  ti0  +  2nN/w)  =  v, (z;  xio,  yi0 ,  /,„)  for  integer  N.  The  discrete  sum  over  particles  can  be 
replaced  by  an  integration  over  initial  conditions,  and  we  may  rewrite  (7)  in  the  form 

8J(jr,  t)  =  —  enb\zo  If  dxQdyqiTL  (*0,  >>o)  f  dt0  <rn{t0) 


.  8[r  t (z ,  x0i  y0,  t0)] 

x  Wz;x0,y0, /0)  |Vz(z;  w<m,o)| 


(9) 


where  vz0  is  the  initial  axial  velocity,  Ag  =  irRg  is  the  cross-sectional  area  of  the  waveguide, 
T  =  L/vzo,  and  o-j.fx,,,  y0)  and  o-n(r0)  describe  the  distribution  of  the  initial  conditions  subject  to  the 
normalization 


f  f  dx0dy0  <rL(x0,  y0)  =  Ab,  (10) 

Ag 

/t n  ,  v 

_r/2  rfr0o-"(r0)  -  r,  01) 

where  Ab  is  the  cross-sectional  area  of  the  electron  beam. 


Substitution  of  the  microscopic  fields  and  the  source  current  density  into  Maxwell’s  equations  for 
the  TE  mode  yields 


where  8a,„  =  ehAjmc 2,  /320  =  v20/c,  u>l  =  4ire2nb/m,  (v,,  v2)  are  the  transverse  components  of  the 
electron  velocity  relative  to  the  basis  vectors  e\  -  ex  cos kwz  +  ey  sin<cwz  and  e7  —  —  ex  sin/cwz 
+  ey  coskwz.  For  the  TM  mode,  we  find  the  similar  result 


•  r. 


~7T  8fl/»  +  1  +  *7 T  ~  k2  —  k2„  8fl|n  -  — f/JjoMn 

dz 2  k 1  c  c 


yv,r/+»  +  v2^/+>  +  2%-  v2Ji(ki„r)  sina,\ 

\ - kT - >4) 


t  +  T  +  T  8fl'"  ’  7^"'" 


^ V|  H''/-’  -  v27'/-)  +  2— p  \2J,(klnr)  cos  a) 


In  the  preceding  equations,  Hln,  7}(±\  and  HV±)  are  mode  (i.e.,  polarization)  dependent  quantities 


defined  as 


— - - — - — t~  ;  TE,„  mode 

(X,„2  -  I1)  Ji2(X'„ ) 

H In  —  1 

7^7 — -  ;  mode 

h  (*ln) 


F/±]  sin0/  +  G/*'  cos  ill,  ;  TE,„  mode 

Ti  —  ptf)  CQS  ^  _  Qi*)  sjn  ^  moc|e 


W,1*'  = 


where 


Fi*'  cosi hi  —  Gi,f>  sin  i /»;  ;  rf(n  mode 
(F/±(  sin  ifr/  +  G;(±l  cos  <bf)  ;  rM/„  mode 


'll/  =  *0  +  X  dz  >  k  +  lk*  ~  V" 

is  the  phase  relative  to  the  ponderomotive  frame,  ib0(=  —  <nt0)  is  the  initial  phase. 


F/*'  =  J,-\(k,„r)  cos  (/ -  l>x  ±  Ji+\(ki„r )  cos  (/  +  I >x - 


Gf(±l  =  J,-\(k,„r)  sin  (/  -  1)*±  J,+X(kinr)  sin  (/  +  1)*. 


Finally, 


<F>  =  —4 


dtocrMo)  S!  d90dr0r0(Tl  (r0,  0O)F 


>v  y  •••  ■■  .v--. . 


Y-Y-V-Y Y-Y-Y-Y  .-•Y-Y-Y-Y-  ,--y v-y-y -y-y  ,-.v.  y-V-YY.-.Y-Y-Y 

\-  -  Vv  ‘I- ‘•’•'•/•Y-Y Y-Y-Y-Y^  V*  -Y-Y-Y-Y- Y-Y- Y-Y-Y-Y-' YY'-Y 

\*  •,*  %"  \  \*  %’  *,  S  V*  «.*  -  *  *■  «.•  »  *  «  '  »•  O  f  m  *  »  •  «  *  »  *  •v'*  .’*«.*  «  •<  i  »>\  V*  •  • 
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describes  Che  average  of  the  beam  electrons  over  both  the  axial  phase  and  cross-section  of  the 
waveguide.  It  is  important  to  recognize  that  this  average  includes  the  effect  of  the  overlap  of  the  elec¬ 
tron  beam  with  the  transverse  mode  structure  of  the  radiation  Held  (often  included  in  one-dimensional 
formulations  in  an  ad  hoc  manner  by  the  inclusion  of  a  filling  factor)  in  a  self-consistent  way. 

In  order  to  complete  the  formulation,  the  electron  orbit  equations  in  the  presence  of  the  static  and 
fluctuation  fields  must  be  specified.  Since  we  describe  an  amplifier  configuration,  we  choose  to 
integrate  in  z  and  write  the  Lorentz  force  equation  in  the  form 

v,-£p  =  -  e8Ei„  -^»x(B0+B,  +  SB,,),  (23) 

where 

S£/,  =  -  ;8B/b  =  V  x  SA/b. 

Substitution  of  the  appropriate  form  for  the  vector  potential  shows  that 

v.-j P\  =  — ^(fl0  -  ykw\z  +  2Clwll(\)  sin  x>/>2  +  -n,ft/2(X)  sin  2x 
uz  y  y 

-  —  mcSa,B[(w  -  kvz)  I V,(~)  -  2 k,n\1J,(k,„r)  cos a,  -  r,„vz  r/+)j,  (24) 

vz-~ p2  -  —  (n0  -  y k„vz  +  2n„/,(\)  sin*)pi  — -fl»pz[/0(\)  +  /2(\)  cos2xl 

dz  y  y 

+  -j  me  5a,fl|(tti  -  fcv2)r,(-)  -  2kln\\J,(k,nr)  cos a,  +  r,„vz  If,(+)j,  (25) 

v- ~-pz  =  —  nH,p2l/0(X)  +  /2(X)  cos2xl  — —  n„p,/2(X)  sin2x 

dz  y  y 

mc8aln\k(vitVl(->  -  v2r,(->)  +  r,„(v,r/+)  +  v2ff',(+')|,  (26) 

for  the  TE,„  mode,  where  n0w  =  \eB0w/mc\,  y  =  (1  —  v2/c2)“1/2  is  the  relativistic  factor,  p  is  the 
momentum,  and 

represents  the  growth  rate  of  the  wave  mode.  For  the  TM,„  mode,  we  find 


(27) 


v, — p i  - (ft0  -  ykwvz  +  2ft  „/i(X)  sinx)p2  +  — ft  h,/»j/2(\)  sin2x 

dz  y  y 


1 


-  -me  5 a,„ 


k2  +  k,2n 


r ,„v2r/ 


y,~~P2  =  —  (fto  “  ykwvz  +  2ft  w/i(X)  sinx>Pi -ft  „p.  l/o(X)  +  /2(X)  cos2x) 

'  dz  y  y 


+  -mcha,H 


k 2  +  *,2 


7-/-’  +  r„v,  if; 


<  +  > 
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In  a  addition,  we  have  that  for  either  the  TE/n  or  TM,n  modes 


v. -^-x  =  v,  cos  kKz  —  v2  sin  kH,z , 


</ 

v2— -.y  =  v,  sin  ac„z  +  v2  cos  kwz , 
dz 


and 
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Both  the  linear  and  nonlinear  evolution  of  the  FEL/Ubitron  amplifier  are  included  in  the  formula¬ 
tion  through  Eqs.  (12MI5)  for  the  fields,  and  (23M32)  for  the  particles. 


III.  NUMERICAL  SIMULATION 


The  set  of  coupled  differential  equation  derived  in  Sec.  II  is  solved  numerically  for  an  amplifier 
configuration  in  which  a  single  wave  of  frequency  w  is  injected  into  the  system  at  z  =  0.  Maxwell's 
Eqs.  ( 1 2 ) - ( 1 5 )  can  be  reduced  to  a  set  of  three  first  order  differential  equations  for  5a, „,  r,„,  and  k. 
Hence,  the  numerical  resolution  of  the  problem  consists  in  the  simultaneous  solution  of  6NT  +  3  first 
order  ordinary  differential  equations,  where  Nr  is  the  total  number  of  electrons.  The  algorithm  we 


employ  involves  the  use  of  a  4th  order  Runge-Kutta  method  to  calculate  the  first  three  steps  after  the 
initial  state,  after  which  an  Adams-Moulton  predictor/corrector  is  employed  for  all  further  steps.  The 
averages  in  Eqs.  (1 2)-(l 5)  are  performed  by  means  by  an  Nth  order  Gaussian  quadrature  technique  in 
each  of  the  variables  r0,  0O);  hence,  N3  =  NT.  For  all  cases  discussed  in  this  work  a  choice  of 
N  =  10  was  found  to  provide  an  accuracy  of  better  than  .1%. 


The  initial  state  was  chosen  to  model  the  injection  of  a  solid,  axi-symmetric  monoenergetic  elec¬ 
tron  beam  of  zero  emittance  and  uniform  cross-section.  Hence,  we  choose  <rL  =  an  =  1,  and  initially 
set  pL  =  0  and  pM  =  mc-Jy $  -  1 .  In  addition,  the  electron  positions  are  chosen  b,  means  of  the  Gaus¬ 
sian  algorithm  within  the  ranges  —it  ^  4i0  <  n,  0  ^  0O  <  2w,  and  Rm{„  <  r0  <  hence,  we  may 

model  the  case  of  either  a  solid  or  an  annular  electron  beam.  Within  the  context  of  this  beam 
geometry,  the  plasma  frequency  is  related  to  the  total  beam  current  !b  by  means  of  the  relation 
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In  order  to  self-consistently  satisfy  the  neglect  of  space-charge  fields  and  treat  the  High-Gain  Compton 
regime,  we  must  require  that 
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within  the  uniform-wiggler  region.  This  limitation  on  the  beam  current  will  vary  with  the  choice  of 


(34) 


energy,  wiggler  and  axial  magnetic  field  strengths,  and  wiggler  period,  but  will  be  adhered  to  in  all  the 
cases  we  study.  The  initial  conditions  on  the  radiation  field  are  chosen  such  that  I'M(r  =  0)  -  0  and 
k{z  =  0)  =  —  k2„  for  an  arbitrary  initial  power  level  Observe  that  the  time-averaged  Poynting 


flux,  Pw,  for  the  waveguide  mode  is  related  to  the  field  amplitude  by  the  relation 


for  the  TE,„  mode,  and 


for  the  TMtn  mode. 
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It  is  important  to  observe  that,  although  the  analysis  applies  to  the  High-Gain  Compton  regime, 
collective  effects  are  included  through  the  dielec',  ic  response  uf  the  plasma  to  the  waveguide  mode, 
and  the  analysis  is  not  purely  single-particle.  Thus,  while  the  wavenumber  is  set  initially  by  means  of 
the  vacuum  dispersion  relation,  the  system  evolves  through  an  initial  transient  regime  into  a  fully  self- 
consistent  dielectrically  loaded  waveguide  mode. 


The  first  case  we  consider  is  that  of  a  wide-band  35  GHz  Ubitron  amplifier  operating  in  the  TEn 
mode.  To  this  end,  we  assume  a  wiggler  field  with  an  amplitude  and  period  of  Bw  —  2kG  and 
X  =  1.175  em,  and  a  waveguide  radius  of  Rg  -  0.36626  cm.  The  electron  beam  configuration  is  that  of 
a  solid  (i.e.,  pencil)  beam  with  an  energy  of  250  keV,  a  current  of  35  Amperes,  and  an  initial  radius  of 

“  0.155  cm.  The  initial  distributions  in  the  axial  phase  space  and  beam  cross-section  are  shown  in 
Figs.  1  and  2.  Each  dot  in  the  illustration  of  the  axial  phase  space  (Fig.  1)  describes  a  "phase-sheet" 
composed  of  100  electrons  distributed  throughout  the  cross-section  of  the  beam.  Each  phase-sheet, 
therefore,  represents  a  cross-sectional  slice  of  the  beam,  which  is  chosen  initially  as  shown  in  Fig.  2. 
The  circle  shown  in  Fig.  2  represents  the  waveguide  wall.  Each  phase-sheet  is  initially  chosen  to  be 
identical;  however,  the  subsequent  evolution  of  the  particle  trajectories  in  the  presence  of  the  radiation 
field  is  followed  self-consistently.  It  should  be  remarked  that  the  particle  distribution  described  herein 
represents  a  uniform  electron  beam.  The  positions  in  (r0,  0O,  0to)  were  chosen  by  a  10-point  Gaussian 
weighting,  and  the  nonuniformity  in  the  positions  of  the  electrons  is  compensated  for  by  a  nonuniform 
weighting  of  the  particles. 

We  now  digress,  briefly,  to  describe  the  types  of  single-particle  trajectory  in  the  combined  wiggler 
and  axial  guide  magnetic  fields.  The  optimum  orbit  for  FEL/Ubitron  operation  is  characterized  by  a 
coherent  electron  motion  in  phase  with  the  wiggler.  Such  orbits  would  have  a  constant  axial  velocity 


Beam  Cross-Section  (kwz  =  0) 


>vv. 


(to  preserve  the  resonance  condition  between  the  particles  and  the  ponderomotive  wave  formed  by  the 
beating  of  the  wiggler  and  radiation  fields)  as  well  as  a  transverse  velocity  of  constant  magnitude.  A 
class  of  orbits  of  this  type  may  be  found  in  the  presence  of  both  a  wiggler  and  an  axial  guide  field. 2S_27 
These  ideal  orbits  describe  helical  trajectories  about  the  axis  of  symmetry  (axi-centered) ,  and  are 
characterized  by 

v  -  vwe!  +  V||e,,  (37) 

X  —  +  v„/vn,  and  x  —  ±  jt/2,  where  (vw,  V||)  are  constants,  and 

v  2n,vn/1(x)/x 

*  n„-  y*wV||  ±  2nw/,(x)  •  ' 

Self-consistent  solution  for  the  axial  velocity  as  a  function  of  B0,  Bw,  and  y  may  be  found  by  sub¬ 
stitution  of  the  expressions  for  y  and  v„  into 

+  v»  -  0  -  y f2-  (39) 

The  solution  is  obtained  by  numerical  means,  and  is  shown  in  Fig.  3  as  a  function  of 
B0  for  Bw  -  2kG,  -  1.175  cm,  and  Vb  —  250  keV  (y  =  1.489).  As  shown  in  the  figure,  there  are 
two  types  of  orbit.  Group  I  orbits  are  found  for  relatively  low  axial  guide  fields  and  are  characterized 
by  ft0  <  ykwv\\.  In  contrast.  Group  II  orbits  occur  for  high  axial  fields  and  have  fl0  >  y*wvi|-  The 
dashed  line  represents  orbitally  unstable  trajectories. 

The  purpose  of  the  adiabatic  entry  taper  in  the  wiggler  field  is  to  inject  electrons  into  the  wiggler 
region  onto  trajectories  approximating  these  ideal  helical  orbits.  Numerical  integration  of  the  particle 
trajectories  (i.e.,  in  the  absence  of  a  radiation  field)  shows  that  such  injection  is  feasible  for  particles 
initially  on-axis  (i.e.,  r(z  =  0)  —  0)  as  long  as  the  orbit  parameters  are  not  too  dose  to  the  magnetic 
resonance  at  fl0~  yX»V||.  Thus,  there  are  practical  limitations  on  the  operation  of  the  Ubitron/FEL 
near  gyroresonance.  Of  course,  the  injection  process  is  not  perfect,  and  the  orbits  differ  slightly  from 
the  ideal  trajectories  in  that  periodic  motion  is  also  observed  corresponding  to  Larmor  oscillations  due 
to  the  axial  guide  field,  Betatron  oscillations  due  to  the  transverse  gradient  in  the  wiggler  field,  and 
higher  harmonics  of  the  wiggler  period.  These  problems  intensify  somewhat  for  electrons  which  are 
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Group  I  Orbits 


injected  off-axis  (i.e.,  r(z  —  0)  >  0),  in  addition  to  which  the  magnitudes  of  the  transverse  and  axial 
velocities  of  the  ideal  helical  component  of  the  motion  vary  with  radial  position.  This  occurs  because 
the  electrons  "see”  a  higher  average  wiggler  field  as  the  displacement  from  the  axis  grows,  and  an  elec¬ 
tron  injected  at  an  initial  radius  of  r0  will  be  characterized  by  a  wiggler  velocity  of28 

2fl  wV||/0(A.o)/  i(X)/X. 

V  sas  — —  -  ■  -  —  ■  ■  — ■  ( 40 1 

*  n0-y*wV||  ±  2fi W/0(X0)/,(X)  ’ 

where  X0  =  fc»,r0.  As  a  consequence,  the  finite  extension  of  the  electron  beam  will  correspond  to  a 
spread  in  (v„,  V||)  across  the  beam  which,  in  turn,  will  result  in  a  broadening  of  the  wave-particle  reso¬ 
nance.  The  net  result  of  all  of  these  effects  is  that  the  beam  may  display  a  very  complex  overall  bulk 
motion,  and  we  emphasize  that  all  of  these  effects  are  included  in  the  simulation  in  a  self-consistent 
way. 

In  order  to  illustrate  the  particle  trajectories  which  occur  in  the  combined  field  structure,  it  is  use¬ 
ful  to  show  the  evolution  of  the  beam  cross-section  at  a  series  of  axial  positions  within  the  interaction 
region.  We  consider  a  set  of  parameters  consistent  with  Group  I  orbits  and  choose  an  axial  field  of 
B0  -  1.3  kG.  This  is  sufficiently  far  from  the  gyroresonance  (see  Fig.  3)  that  injection  onto  near-ideal 
orbits  is  possible.  In  Fig.  4,  we  show  the  beam  cross-section  at  kwz  =  75,  which  corresponds  to  a  point 
just  after  the  start  of  the  uniform  wiggler  region  at  kwz  =  62.83.  Note  that  the  initial  value  of  the  radi¬ 
ation  field  was  chosen  to  correspond  to  an  input  power  of  P,„  —  10  W  at  a  frequency  ut/ckw  -  1.3,  so 
that  the  perturbation  to  the  single-particle  orbits  at  this  point  due  to  the  electromagnetic  field  is  negligi¬ 
ble.  It  is  immediately  evident  from  the  figure  that  the  beam  has  been  substantially  compressed  due  to 
the  focussing  effect  of  the  wiggler.  In  addition  the  beam  center  is  shifted  off-axis  corresponding  to  the 
helical  motion  imposed  by  the  wiggler.  Thus,  the  beam  has  "spun-up”  due  to  the  transverse  velocity 
arising  from  the  combined  influence  of  the  wiggler  and  guide  magnetic  fields.  The  injection  process  for 
these  parameters  yields  orbits  close  to  the  ideal  helical  trajectories,  and  kwrceme,  —  .25  for  the  beam 
center  which  is  in  good  agreement  with  the  result  expected  from  the  ideal  trajectories.  The  beam  also 
displays  a  rotational  motion  about  the  center  which  is  accounted  for  by  the  combined  effects  of  the 


variation  in  vw  across  the  beam  cross-section  and  the  betatron  oscillation  due  to  the  transverse  wiggler 
gradient. 

Overall  bulk  motion  of  the  beam  is  expected  to  display  a  periodicity  at  the  fundamental  wiggler 
harmonic,  and  to  twist  the  beam  into  a  helix  about  the  axis  of  symmetry.  Such  motion  is  observed  as 
is  shown  in  Figs.  5-8,  which  displays  the  beam  cross-section  as  it  evolves  from  kwz  “  150  through 
kwz  —  156  (approximately  one  wiggler  period).  Note  that  this  is  well  within  the  region  of  strong  linear 
(i.e.,  exponential)  growth  of  the  radiation  field,  and  substantial  amplification  of  the  input  signal  has 
occurred.  As  a  result,  each  phase-sheet  has  reacted  in  a  slightly  different  way.  This  is  indicated  by  the 
"smearing"  of  the  dots  in  Figs.  5-8  which  marks  small  displacements  of  the  beam  centers  in  each  phase 
sheet  as  well  as  small  differences  in  the  rates  of  rotation  about  the  beam  centers. 

We  conclude,  therefore,  that  the  motion  of  the  beam  includes  an  overall  bulk  helical  motion 
about  the  symmetry  axis  at  the  wiggler  period,  as  well  as  a  "pinwheeling"  type  of  rotation  about  the 
beam  center  at  a  much  longer  period. 

The  evolution  of  the  waveguide  mode  corresponding  to  these  transverse  beam  motions  is  shown 
in  Fig.  9.  The  electromagnetic  wave  was  a  TEU  waveguide  mode  (for  Rg  —  0.36626  cm)  at 
m/ckw  «*  1.3  which  corresponds  to  a  frequency  of  /  =  33.2 GHz.  Beam  voltage  and  current  were  250 
kV  and  35  A  respectively  for  an  average  beam  power  of  8.75  MW.  As  shown  in  the  figure,  the  growth 
of  the  wave  mode  was  approximately  exponential  after  an  initial  transient  period  for  kwz  <  80.  During 
the  linear  phase  of  the  interaction  the  growth  rate  was  T\\fkw  =  0.029,  and  a  small  increase  was 
observed  prior  to  saturation  at  kwz  -  267.  The  radiation  power  at  saturation  was  1.87  MW  for  an 
overall  efficiency  of  rj  =  21.4%.  A  complete  spectrum  of  the  efficiency  at  saturation  versus  frequency 
is  shown  in  Fig.  10,  in  which  the  dots  represent  the  numerical  results  of  the  simulation.  We  observe 
that  peak  efficiency  occurs  for  ~  1.3c/cw,  with  a  bandwidth  Aw/a> ^  ~  54%. 
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Saturation  is  by  means  of  particle  trapping  in  the  ponderomotive  wave  formed  by  the  beating  of 
the  radiation  and  wiggler  Helds.  This  is  clearly  shown  in  Pig.  1 1  in  which  we  plot  the  axial  phase  space 
at  saturation  for  the  case  of  <u/ckw  -  1.3.  The  solid  lines  in  the  figure  represent  approximate  separa- 
trices  calculated  under  the  assumption  that  all  particles  are  executing  the  ideal  helical  trajectories. 
Thus,  the  actual  number  of  electrons  trapped  may  differ  slightly  from  that  shown  in  the  figure,  but  the 
conclusion  remains  valid.  The  beam  cross-section  at  saturation  is  shown  in  Fig.  12.  Although  the  fig¬ 
ure  seems  to  show  a  chaotic  state  in  which  the  initial  uniformity  in  cross-section  of  each  phase-sheet 
has  been  destroyed,  a  closer  examination  of  the  positions  of  the  electrons  within  each  phase-sheet 
shows  that  the  electron  positions  have  retained  their  coherence  and  appear  much  as  shown  in  Figs.  5-8. 
However,  the  beam  centers  in  each  phase-sheet  have  shifted  markedly  so  that  the  superposition  of  the 
electrons  in  all  the  phase-sheets  gives  the  appearance  of  a  random  distribution. 

The  consistency  of  the  code  has  been  checked  several  ways.  The  most  fundamental  is  the 
requirement  of  energy  conservation  between  the  particles  and  the  wave,  and  agreement  between  the 
energy  lost  by  the  beam  and  that  gained  by  the  wave  was  found  to  be  significantly  better  than  0.1%  in 
all  cases  considered.  In  addition,  the  growth  rate  in  the  linear  regime  has  been  checked  against  the 
predictions  of  a  three-dimensional  linear  theory  of  the  interaction.  Although  the  linear  theory 
represents  an  idealized  model  in  which  all  particle  trajectories  are  described  by  the  ideal  helical  orbits, 
agreement  between  the  simulation  and  the  linear  theory  is  good.  At  the  frequency  corresponding  to 
peak  efficiency  Wcfc*  —  1.3),  the  linear  theory  predicts  a  growth  rate  of  r^n/Ac^  =  0.031,  which  com¬ 
pares  well  with  the  simulation  result  of  Tu/k*  =  0.029.  Finally,  an  estimate  of  the  efficiency  on  the 
basis  of  the  simple  heuristic  phase-trapping  model  yields  an  estimate  of  a  19.1%  efficiency,  which  is 
close  to  the  21.4%  efficiency  obtained  from  the  simulation. 

The  case  of  the  TMU  mode  has  also  been  examined  using  the  simulation  code.  In  order  to  facili¬ 
tate  comparison  with  the  TEU  mode,  the  parameters  specifying  the  external  magnetic  fields  and  elec¬ 
tron  beam  remain  unchanged;  that  is  B0—\.3kG<  Bw~2.0kG,  A w  —  1.175  cm ,  Vb  —  250  kV, 
4-35/4,  and  the  initial  beam  radius  is  Rb  -  0.155  cm.  The  only  alteration  is  in  the  choice  of  the 
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waveguide  radius  in  order  to  shift  the  waveguide  dispersion  curve  relative  to  the  beam  resonance  line 
(i.e.,  o>  —  (k  +  kw)v\\).  We  choose  Rg  -  0.76223cm  so  that  the  waveguide  cutoff  is  identical  to  that  of 
the  TEU  mode,  and  the  interaction  occurs  over  the  same  frequency  band.  Numerical  results  then  yield 
a  peak  efficiency  T)max  =  5.92%  at  a  frequency  of  to/ckw  =  1.78.  This  is  significantly  lower  than  the 
peak  efficiency  of  21.4%  found  for  the  TEU  mode  at  a  frequency  of  u/ckw  =  1.3.  The  growth  rate  was 
also  significantly  lower,  and  we  observe  a  growth  rate  of  only  r \\/kw  =  0.012  at  the  frequency 
corresponding  to  peak  efficiency.  For  this  reason,  we  shall  confine  our  attention  in  the  remainder  of 
this  work  to  the  TE  modes.  However,  we  emphasize  that  no  attempt  was  made  to  optimize  the  interac¬ 
tion  with  respect  to  the  TM\\  mode,  and  it  would  be  unjustified  to  conclude  (as  in  the  case  of  the 
Gyrotron,  for  example)  that  the  FEL/Ubitron  interaction  favors  the  TE  mode. 

The  cases  discussed  thus  far  dealt  with  waveguide  modes  with  an  /  —  1  azimuthal  mode  number. 
Examination  of  the  dynamical  equations  shows  that  as  long  as  \  is  an  approximate  constant  along  the 
particle  trajectories  (\  —  ±  ir/2  for  the  ideal  helical  orbits)  a  selection  rule  exists  whereby  TEh  or  TM,„ 
modes  interact  via  wave-particle  resonances  at  w  =  (k  +  lkw)v\\.  It  has  been  conjectured,  therefore, 
that  no  interaction  is  possible  for  the  /  -  0  modes,  and  numerical  study  of  test  cases  for  TE0i  and  TM0l 
modes  has,  indeed,  shown  that  no  gain  occurs.  Higher  harmonic  interactions  may  also  be  studied  by 
the  selection  of  higher  order  (i.e.,  /  >  1)  modes,  and  the  results  of  a  series  of  simulations  is  shown  in 
Fig.  13  for  the  TE2 1  mode.  As  in  the  case  of  the  TMU  mode,  the  parameters  specifying  the  external 
magnetic  fields  and  the  electron  beam  are  chosen  to  be  identical  to  those  used  to  study  the  TEn  mode 
in  order  to  facilitate  comparison  of  the  interaction  at  the  two  harmonics.  Only  the  waveguide  radius  has 
been  changed  to  R„  -  0.28  cm  in  order  to  bring  the  intersection  points  between  the  TE2\  dispersion 
curve  and  the  beam  resonance  line,  &>«(*  +  2*„)v||,  sufficiently  close  together  to  obtain  a  relatively 
broad-band  gain  spectrum.  The  growth  rate  found  at  the  frequency  which  corresponds  to  peak  effi¬ 
ciency  (i.e.,  o»p KtJckw  =  2.9)  was  r2)/fcw  =  0.025  which  is  comparable  to  that  found  for  the  TEU  mode 
(fu/^»  —  0.029  at  u>/ckw  =  1.3).  The  efficiency  is  plotted  as  a  function  of  frequency  for  the  TE2\ 
mode  in  Fig.  13.  While  the  bandwidth  Aw/w,*.*  —  43%  is  comparable  and  the  frequency  is  more  than 
doubled  relative  to  the  TE\\  mode,  the  peak  efficiency  of  4.7%  is  greatly  decreased. 
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We  now  turn  to  the  case  of  a  higher  energy  electron 'beam.  In  Fig.  14  we  plot  the  variation  of  the 
axial  velocity  of  the  ideal  helical  orbits  versus  B0  for  1  MeV  electrons  and  a  wiggler  characterized  by 
Bw  —  1.0  AG  and  A*  —  3  cm.  As  indicated  in  the  figure,  Group  I  orbits  are  found  for  axial  guide  fields 
Bq  <  6 ,4kG,  and  Group  II  orbits  for  B0  >  8.2 kG.  No  such  orbits  are  possible  between  these  limiting 
values.  The  first  case  we  consider  is  that  of  B0  —  3.5  which  corresponds  to  a  Group  I  type  of  trajectory, 
and  we  assume  a  1  MeV,  50  Amp  electron  beam  with  an  initial  radius  of  0.2  cm.  The  waveguide  radius 
was  chosen  to  be  Rg  =  0.45  cm,  so  the  upper  and  lower  intersection  frequencies  between  the  TEU 
dispersion  curve  and  the  beam  resonance  line  are  well-separated  and  occur  at  ai/ckw  —  10.9  and  2.8 
respectively.  Gain  is,  indeed,  found  in  the  vicinity  of  these  frequencies,  as  shown  in  Fig.  15  in  which 
we  plot  the  interaction  efficiency  versus  frequency  for  this  case.  We  observe  that  the  efficiency  drops 
rapidly  near  the  upper  (lower)  range  of  the  lower  (upper)  frequency  band.  The  lower  frequency  band 
corresponds  to  frequencies  /  —  28  GHz  which  are  comparable  to  those  found  for  the  TEU  mode  for  a 
250  keV  beam,  and  although  the  bandwidth  is  narrower  than  that  found  previously  (see  Fig.  10)  the 
efficiencies  are  comparable.  The  upper  frequency  range  shows  a  broader  bandwidth  and  somewhat 
reduced  efficiency  (tj,*^  —  13.7%). 

The  results  of  a  series  of  simulations  for  B0  —  1 1.75  AG  and  Rg  —  0.5  cm  are  shown  in  Fig.  16  for 
parameters  consistent  with  Group  II  orbits,  in  which  we  plot  the  interaction  efficiency  as  a  function  of 
frequency.  It  should  be  remarked  that  for  some  of  these  runs  a  small  fraction  of  the  beam  was  lost  to 
the  wall.  The  procedure  followed  when  this  occurred  was  to  eject  such  particles  from  the  simulation, 
and  the  efficiencies  shown  include  this  effect.  The  bandwidth  observed  was  Aw/«Opeak  —  57%  over  a 
frequency  range  of  liGHz  -  45 GHz.  However,  while  the  frequency  range  is  comparable  to  that  used 
for  the  lower  energy  (250  keV)  case,  the  observed  efficiencies  are  considerably  higher.  We  observe  a 
peak  efficiency  rtpeak  —  47%  at  frequencies  in  the  range  =  3.0-3. 1,  which  is  more  than  dou¬ 

ble  that  found  for  the  Group  I  orbits  with  a  1  MeV  beam  in  this  frequency  range.  It  should  be 
emphasized  here  that  this  efficiency  is  found  with  uniform  wiggler  and  guide  magnetic  fields,  and  no 
magnetic  field  tapering  was  used.  This  enhancement  in  the  interaction  efficiency  associated  with  Group 
II  orbit  parameters  has  also  been  observed  in  one-dimensional  simulation. w 
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Fig.  15  —  Spectrum  of  the  interaction  efficiency  versus  frequency  for  the  TEU 
mode  and  parameters  corresponding  to  Group  I  orbits 
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The  small  degree  of  variation  in  the  efficiency  about  the  line  of  best  fit  is  difficult  to  precisely 
account  for  due  to  the  complexity  of  the  system  of  coupled  equations.  Some  of  the  causes  of  this  addi¬ 
tional  structure  undoubtedly  include:  (1)  variations  in  the  number  of  particles  which  are  lost  to  the 
wall,  and  the  positions  (in  z )  at  which  they  are  lost,  (2)  the  variation  in  the  ponderomotive  potential 
with  the  frequency  at  which  the  interaction  occurs,  and  (3)  variations  in  the  evolution  of  the  radial  pro¬ 
file  of  the  beam  near  saturation  (where  the  radiation  field  is  large)  have  an  effect  upon  both  the  overall 
filling  factor  and  the  trapping  fraction.  However,  this  is  an  area  of  ongoing  investigation,  and  we  can¬ 
not  rule  out  the  possibility  of  an  additional  modulation  of  the  interaction  due  to  departures  of  the  elec¬ 
tron  trajectories  from  the  ideal  helical  orbits. 

IV.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  developed  a  fully  self-consistent  nonlinear  theory  and  numerical  simulation 
of  the  Ubitron/FEL  amplifier  in  three-dimensions.  The  particular  configuration  considered  consists  of 
a  cylindrical^  symmetric  electron  beam  of  arbitrary  cross-section  injected  into  a  loss-free  cylindrical 
waveguide  in  the  presence  of  both  a  helically  symmetric  wiggler  and  a  uniform  axial  guide  magnetic 
field.  In  addition,  the  adiabatic  injection  of  the  electron  beam  has  been  modeled  by  including  an  initial 
taper  of  the  wiggler  field  amplitude.  The  system  of  equations  derived  describes  the  self-consistent  evo¬ 
lution  of  both  the  wave  fields  and  the  trajectories  of  an  ensemble  of  electrons.  The  analysis  has  been  • 

performed  for  both  the  TE  and  TM  modes,  and  describes  the  overlap  of  the  transverse  mode  structure 
and  the  electron  beam  (i.e.,  the  filling  factor)  in  a  self-consistent  manner.  Space-charge  fields  have 
been  neglected  in  the  analysis,  so  the  treatment  is  applicable  to  the  High-Gain  Compton  regime  of 
operation.  In  addition,  self-field  effects  of  the  electron  beam  have  been  neglected.  Since  the  problem 
of  interest  in  the  Ubitron/FEL  amplifier,  only  single  wave-mode  propagation  is  considered.  This  per¬ 
mits  an  average  over  the  wave  period  to  be  performed  which  eliminates  the  fast  time  scale  phenomena 
from  the  formulation,  and  results  in  a  great  increase  in  computational  efficiency  over  a  full-scale 
particle-in-cell  simulation  code. 


33 


The  electron  trajectories  are  integrated  using  the  exact  Lorentz  force  equations,  so  that  we  are 
able  to  examine  the  detailed  motion  of  the  electron  beam  in  the  combined  field  structure.  Overall  bulk 
motion  of  the  beam  exhibits  a  dominant  oscillation  at  the  fundamental  wiggler  period,  as  expected, 
which  twists  the  beam  into  a  helix  about  the  axis  of  symmetry.  In  addition,  a  significant  focussing 
effect  is  observed  to  occur  due  to  the  radial  inhomogeneity  of  the  wiggler  which  can  result  in  a  signifi¬ 
cant  contraction  in  the  beam  radius.  Superimposed  on  this  bulk  motion  is  (1)  a  shear  in  the  electron 
velocity  across  the  beam,  and  (2)  a  slow  time-scale  Betatron  oscillation  whose  net  effect  results  in  a 
slow  pinwheeling  motion  about  the  beam  center.  Larmor  motion  due  to  the  axial  guide  field  is  more 
difficult  to  identify  because  the  period  may  be  comparable  (especially  near  gyroresonance)  to  the 
wiggler  period.  However,  the  adiabatic  beam  injection  technique  was  found  to  be  effective  as  long  as 
the  parameters  were  not  too  close  to  gyroresonance,  and  the  Larmor  motion  can  be  kept  small. 

Numerical  simulations  were  conducted  for  parameters  corresponding  to  a  250  keV  Ubitron  and  a  1 
MeV  FEL.  In  all  cases,  comparison  of  the  energy  lost  by  the  particles  to  that  gained  by  the  wave 
showed  energy  to  be  conserved  to  within  an  accuracy  significantly  better  than  0.1%.  The  consistency  of 
the  simulation  was  also  checked  by  comparison  of  (1)  the  growth  rate  found  by  the  code  with  that 
predicted  from  a  linear  theory  of  the  instability,24  and  (2)  the  saturation  efficiency  obtained  in  the 
simulation  with  that  found  by  simple  phase-trapping  arguments.  In  both  cases,  good  agreement  was 
found. 

The  parameters  chosen  to  model  the  250  keV  Ubitron  and  the  1  MeV  FEL  were  chosen  so  that 
the  output  frequency  was  in  the  neighborhood  of  35  GHz  for  both  cases.  In  each  case,  the  peak 
efficiency  for  the  TEU  mode  was  found  to  be  in  the  neighborhood  of  20%  for  parameters  corresponding 
to  Group  I  trajectories  (i.e.,  a  relatively  low  axial  guide  field).  However,  the  Group  II  type  of  trajectory 
was  found  to  result  in  substantial  enhancements  in  both  the  efficiency  (—  45%)  and  bandwidth  for  the 
case  of  the  higher  energy  (1  MeV)  beam.  This  is  consistent  with  the  results  found  for  a  previous  one¬ 
dimensional  simulation  of  a  1.25  MeV  beam.  The  amplification  of  the  TMU  was  also  studied  for 
parameters  associated  with  the  250  keV  Ubitron,  and  found  to  result  in  substantially  lower  growth  rates 


and  saturation  efficiencies  than  the  TEU  mode.  However,  since  no  optimization  of  the  parameters  for 
the  TM  mode  was  attempted,  it  is  premature  to  conclude  that  the  FEL/Ubitron  interaction  favors  the 
TE  mode. 

Higher  harmonic  interactions  have  also  been  studied,  and  arise  from  the  azimuthal  variation  of 
the  modes  in  a  cylindrical  waveguide.  Since  the  phase  of  the  TEln  and  TM,„  modes  vary  as 
exp(</fcz  +  ilO  -  (at)  and  the  azimuthal  variation  of  the  ideal  helical  orbits  is  given  by  0  —  kKz ,  the 
beam  resonance  condition  is  <o  —  (/c  +  lkw)v.  and  both  the  TE,n  and  TM,„  modes  resonate  at  the  Ith 
Doppler  upshift.  It  should  be  remarked  that  such  higher  harmonic  interactions  do  not  depend  upon  a 
corresponding  higher  harmonic  component  in  the  single- particle  orbits.  In  view  of  this,  simulation  runs 
were  made  for  the  TE 21  mode,  and  growth  rates  were  found  to  be  comparable  to  those  for  the  TE u 
mode  at  approximately  half  the  frequency.  The  saturation  efficiency,  however,  was  found  to  be  greatly 
reduced  relative  to  the  TEn  mode. 

In  view  of  the  high  efficiencies  found  for  the  35  GHz  examples  shown,  important  future  areas  of 
investigation  include  the  detailed  scaling  of  the  saturation  efficiency  at  a  given  frequency  with  such 
parameters  as  the  beam  energy  and  the  wiggler  parameters.  The  extremely  high  efficiencies  found  for 
parameters  associated  with  high  axial  field  Group  II  trajectories  lend  particular  importance  to  a  determi¬ 
nation  of  the  operational  limits  of  this  regime.  Finally,  it  should  atso  be  noted  that  these  results  have 
been  obtained  for  a  monoenergetic  electron  beam,  and  the  inclusion  of  a  finite  energy  spread  in  the 
simulation  is  currently  being  pursued. 
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